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Abstract 

Two-dimensional CP^" 1 models are used to compare the behavior of different 
cooling techniques on the lattice. Cooling is one of the most frequently used tools to 
study on the lattice the topological properties of the vacuum of a field theory. We 
show that different cooling methods behave in an equivalent way. To see this we apply 
the cooling methods on classical instantonic configurations and on configurations of 
the thermal equilibrium ensemble. We also calculate the topological susceptibility 
by using the cooling technique. 
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1 Introduction 



Two-dimensional (2d) CP^ -1 models JT|, ^ play an important role in quantum field the- 
ory since they represent a convenient theoretical laboratory to investigate analytical and 
numerical methods to be eventually exported to QCD. This is believed to be possible 
since these models have several important properties in common with QCD, for instance 
asymptotic freedom, spontaneous mass generation, confinement and topological structure 
of the vacuum ||. In particular, the CP^ -1 vacuum, like the QCD one, admits instan- 
ton classical solutions. The non-trivial structure of the vacuum for the CP W_1 models, 
as well as for QCD, can be probed by defining suitable operators, like the topological 
charge density and its correlators, and by applying the prescriptions of quantum field 
theory. Yet most of the observables relevant for the investigation of the vacuum topolog- 
ical structure of a theory are beyond the reach of perturbation theory. Therefore some 
approximation is needed which works also in the non-perturbative regime. So far, the 
most effective and comprehensive tool is represented by Monte Carlo simulations on a 
space-time lattice. The lattice approach for the study of the vacuum structure brings 
about a very delicate task, that of separating the quantum fluctuations at the level of the 
lattice spacing from the relevant long-distance topological ones in any given configuration 
of the thermal equilibrium ensemble. Indeed, trying to assign a topological charge to a 
thermal equilibrium configuration, treating on the same footing short-distance (quantum) 
fluctuations and long-distance (topological) ones, would even prevent the reaching of the 
continuum limit. One powerful technique, frequently adopted in the literature to study 
the topological structure of the vacuum, is the so-called cooling method ||. The idea 
behind the cooling method is that the quantum fluctuations at the level of the lattice 
spacing of a given configuration can be smoothed out by a sequence of local minimiza- 
tions of the lattice action, thus leading to a configuration where only the long-distance, 
topologically relevant fluctuations survive. When this procedure is performed on a set 
of well-decorrelated equilibrium configurations, it yields a "cooled" equilibrium ensemble 
on which the expectation values of topological quantities can be determined. The precise 
definition of a cooling procedure is however arbitrary to a large extent: the lattice action 
to be locally minimized can be chosen differently from the one used in the thermalization, 
the number of cooling steps is a free variable, driving or control parameters of the cooling 
can be introduced, and so on. Comparing the behavior of different cooling methods is an 
interesting issue, especially when a new cooling strategy is proposed. 

The aim of the present work is to get insight into a new cooling method first adopted 
in Ref. || for the SU(2) gauge theory and to compare it with two other methods, al- 
ready known in the literature, namely the "standard" Q cooling and its controlled version 
adopted by the Pisa group (see for instance Ref. 0). The new cooling method, as it will 
be pointed out in the next Section, follows a different strategy with respect to the other 
two and could possess various interesting features, inaccessible to the other two cooling 
methods, like for instance preserving instanton-anti-instanton (I-A) pairs with interaction 
energy smaller than a fixable value |5| . The test-field for this comparison will be provided 
by the CP^ -1 models, which have all the interesting topological properties of the SU(iV) 
gauge theories, but are much easier to simulate on a lattice. 
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The paper is organized as follows: in Section we introduce the lattice action and 
the lattice discretization of the relevant topological quantities. In Section [3| we describe 
in detail the three cooling methods. In Sections ^ and [| we discuss their performance on 
1-instanton and I- A classical lattice configurations respectively. In Section [!] we consider 
their behavior on thermal equilibrium configurations and find a correspondence between 
the three cooling techniques. In Section [7] we determine the physical value of the topological 
susceptibility of CP 9 using the new cooling method and compare the results to those from 
an alternative approach not based on cooling. Finally in Section |8] we draw our conclusions. 



2 Lattice definition of action and observables 

The CP^ -1 models describe the physics of a gauge invariant theory of interacting classical 
spins. The continuum action for the 2d CP^ -1 model is 



S = - d 2 x D^z ■ D^z , D ll = d tt + iA ll , (1) 

where g is the coupling constant and z(x) is a A— component complex scalar field which 
obeys the constraint z(x) ■ z(x) = 1. The bar over a complex quantity means complex 
conjugate and a central dot among vectors implies the scalar product: z(x) ■ z(x) = 

We regularize the theory on a square lattice using the standard discretization: 

S L = -N/3 fa + £) • z(x) A M (x) + z(x) ■ z(x + ju) X^x) - 2) , (2) 

x i*=l,2 

where x indicates a generic site of the lattice; A M (x) is a U(l) gauge field (X^x) A At (x) = 1) 
and (3 = l/(Ng) where g is the bare lattice coupling constant. We used the standard 
action both to generate thermal equilibrium configurations in Monte Carlo simulations 
and during the cooling procedure. 

The continuum topological charge is defined as 

Q = J d 2 xq(x) = ^- J (fxe^D^z ■ D v z = ^- j) dx^A^ , (3) 

where e Mi , is the antisymmetric tensor with 612 = 1 and q(x) is the topological charge den- 
sity which can be written as the divergence of a topological current A M (x) = 1/ (271)6^ A u , 
q(x) = <9 M A M (x) |], 0. The last integral in Eq. (Q) is calculated on the 2d plane along the 
circle of infinite radius. 

The topological susceptibility x is a renormalization group invariant quantity which 
gives a measure of the amount of topological excitations in the vacuum. It is defined as 
the two-point zero-momentum correlation of the topological charge density operator q(x), 

X = J d 2 x(0\T[q(x)q(0)]\0) . (4) 
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This equation needs a prescription for the contact term coming from the product of oper- 
ators at the same point. We use the prescription || 

X = j d 2 xd,(0\T[K^x)q(0)}\0) , (5) 

which corresponds to the physical requirement that \ = in the sector with trivial topol- 
ogy (see for instance Ref. |§). 

On the lattice any discretization of the topological charge density with the correct 
naive continuum limit can be used. A standard definition is 

q L (x) = ~ [P(x)A^P(x)A v P(x)} , (6) 

fiu 

with ^ ^ 

^ _ P( x+fi ) - P(x-fj) ^ p ^ x) ^ -, {x)zj{x) _ (7) 



The lattice topological susceptibility is defined as 

X 



i / E^W(o)) = ^((g L ) : 



where Q L ee ^2 x q L (x) and L is the lattice size in lattice units. Unless otherwise stated, 
hereafter the brackets (•) mean average over decorrelated equilibrium configurations. 

In the CP^ -1 model q(x) is a renormalization group invariant operator. Imposing 
that q(x) does not renormalize in the continuum generally implies a finite multiplicative 
renormalization Z on the lattice [jTIJ 

q L {x) = a 2 Z{f3)q{x) + 0(a 4 ) , (9) 

a being the lattice spacing. 

Moreover, the lattice topological susceptibility x L i n general does not meet the con- 
tinuum prescription given in Eq. (|J), thus leading to an additive renormalization which 
consists in mixings of \ L with operators which have the same quantum numbers as Xi 

X L (P)=a 2 Z(P) 2 X + M((3), (10) 

where M{j3) indicates such mixings. The extraction of x from the lattice thus requires the 
determination of both the multiplicative and additive renormalizations, Z(/3) and M(j3). 
This is a feasible task, as will be shown in Section |7|, and in this way one obtains the 
so-called "field theoretical" determination of the topological susceptibility. 

Alternatively one can use cooling to determine the topological susceptibility. This 
method will be described extensively in the next Section. The idea is that, by an iterative 
process of local minimization of the action, the quantum fluctuations at the scale of the 
ultraviolet cut-off, which are responsible for the renormalizations, are suppressed thus 
implying Z — > 1 and M — > in Eq. flT0| ) as the cooling iteration goes on, so that a 2 x 
can be extracted directly from the measurement of x L after cooling, as it will be shown in 
Section [j]. 
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3 The cooling method 



A cooling step applied on a lattice configuration consists in assigning to every field vari- 
able a new value chosen in order to locally minimize the lattice action. The iteration 
of this procedure converges to a configuration which should represent a solution of the 
Euclidean classical equations of motion. In the continuum 2d CP^ -1 models, as well as 
in the case of 4d SU(iV) gauge theories, there are classical solutions that belong to a 
definite topological charge sector, identified by an integer value of the topological charge 
Q and by a finite action (S = 2n\Q\/g in the case of CP^ -1 ). The classical solutions 
with Q 7^ are called instantons or anti-instantons (depending on the sign of Q). On 
the lattice, owing to the artifacts introduced by the discretization of the theory, configu- 
rations with non-trivial topology can represent only approximated solutions of the lattice 
equations of motion. The common belief is that a moderately long cooling sequence can 
be able to wash out the quantum fluctuations of a given lattice configuration, so revealing 
the underlying metastable state with possible non-trivial topology, and that this state is 
the lattice counterpart of a continuum classical one. Of course, if the cooling is further 
protracted, the metastable state will eventually "decay" into the trivial vacuum. Stated 
in other words, once the quantum fluctuations have been eliminated and the topological 
bumps of a configuration (instantons or anti-instantons) have been revealed, if the cooling 
is further protracted then the isolated instantons or anti-instantons shrink, the I-A pairs 
annihilate and the configuration falls into the trivial vacuum of the lattice theory. In 
consideration of the above arguments, it is clear that a great care is needed in determining 
when a cooling iteration has to be stopped. Nevertheless, there is no way to prevent the 
loss of topological signal at lattice scales of the order of the lattice spacing, since topology 
at these scales cannot be distinguished from quantum fluctuations. It must be pointed out 
however, that any loss of topological signal due to cooling which occurs at fixed scales in 
lattice units becomes irrelevant in physical units as far as the coupling is tuned towards its 
critical value, i.e. in the continuum limit, when the lattice spacing a tends to zero. This 
is not the case only when the instanton size distribution is singular in the short-distance 
regime, as in the 2d CP 1 model JXTL [12[. 



In the case of the 2d CP^ -1 models, the cooling algorithm amounts to replacing at 
each lattice site x the field variables z(x) and A M (x) by new ones z new (x) and A^ cw (x), 
chosen in order to locally minimize the lattice action, while the other field variables are 
left unchanged. A cooling step consists in a sweep through the entire lattice volume in 
order to apply this local minimization sequentially at every site x. The terms in the action 
density which depend on z(x) and A M (x) for a given lattice site x have the form 

s z = -2(3NRe{z(x)-F z (x)} , F z (x) = £ A„(z-£) + z{x + $) A„(z)) , (11) 

/U=l,2 

and 

s x = -2pN Re{A M (x) F x {x, //)} , F x {x, n) = z(x) ■ z(x + ju) , (12) 

where F z (x) and F\(x,fi) are the respective forces. It is easy to see that the local min- 
imum of the lattice action is obtained when the new field variables coincide with the 



4 



Figure 1: 6 z (x) (left) and 6\(x, 1) (right) for an "artificial" 1-instanton with size p = 10a 
placed in the middle of a 60 2 lattice. 



corresponding normalized forces: 



r new/ \ _ F z( X ) \new( \_ F \{ X , W // - 1 9 



with = \/u ■ u for a generic complex number or iV-component vector u (z(x), F z (x), 
A Al (x), F\(x,fi), etc.). According to the supplementary conditions under which the above 
replacements are performed, we have different cooling procedures. The first cooling pro- 
cedure, which we call "standard", is the one where the replacements z(x) — > z new (x) 
and A M (x) — > A^ ew (x) are unconstrained. The second one, which we call "Pisa cool- 
ing", is the one where the replacement of z(x) and A M (x) is performed if the distance 
between the old and the new fields is smaller than a previously fixed parameter 5p- lsa : 
\\^ni x ) ~ ^ ew ( x )\ \ — <^Pisa and analogously for the z fields, otherwise the new field is cho- 
sen to lie on the (z(x), z new (x)) ((A M (x), A° ew (x))) plane, at a distance from z(x) (A M (x)) 
exactly equal to 5p isa . The third cooling procedure, called "new cooling" in the follow- 
ing, was introduced in Ref. |J. It consists in making the replacement only if the distance 
between the old and new fields is larger than a fixed parameter 5, otherwise the field is 
left unchanged. We observe that the three methods follow a different strategy. In par- 
ticular, while the Pisa cooling acts first on the smoother fluctuations and the number of 
cooling iterations to be performed has to be fixed a priori, the new cooling performs local 
minimizations only if these fluctuations are larger than a given threshold and the iteration 
stops automatically when all the surviving fluctuations are smoother than the (lattice) 
threshold determined by the dimensionless parameter 5. It is easy to see that the new 
field variables differ from the old ones by terms which, in the continuum limit, are of order 
a 2 , where a is the lattice spacing. Therefore the parameter 5 is expected to scale in the 
continuum limit like a quantity of dimension two. In the following we will compare these 
three cooling methods in several contexts, namely on "artificial" 1-instanton configura- 
tions, on I-A pairs, on thermal equilibrium configurations and in the determination of the 
continuum topological susceptibility. Moreover we will establish a relation between the 
number of cooling iterations in the standard and Pisa coolings with the parameter 5 of the 
new cooling and will discuss the validity of the picture of cooling as a diffusion process. 
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4 Cooling 1— instanton configurations 



In the continuum CP^" 1 models, in the infinite 2d space-time, a 1-instanton configuration 
is a classical solution of the equations of motion characterized by a topological charge one 
and a finite action (S = 2ir/g). The explicit form of a continuum 1-instanton configuration 
is 

w(x 



with 



z ( x ) = I. 7 mi > A* = iz ( x ) ' 9^z(x) , (14) 

umx) = m H f, (15) 

P 

where u and t> are complex iV-vectors which satisfy u • u = v • v — 1 and w • v = 0. 
The space-time coordinates (ci, c 2 ) are the center of the instanton and the real parameter 
p represents a measure of its size. An "artificial" 1-instanton lattice configuration can 
be represented by the discretization of the above continuum configuration on a 2d finite 
volume lattice with periodic boundary conditions. In particular, we used the following 
expressions for the fields z and A J^j 

p x\ — c\- i(x 2 - c 2 ) 



a/p 2 + (si - ci) 2 + (x 2 - c 2 ) 2 ' " a/p 2 + (xi - ci) 2 + (x 2 - c 2 ) 2 



A M (x) 



0, i = 3, ...,iV 

z(x + p) • z(x) 
||^(x + p) • 2(x)|| 



(16) 



which correspond to choose u = (1, 0, 0, . . .) and v = (0, 1, 0, 0, . . .). Actually 1-instanton 
configurations cannot exist on a 2d torus even in the continuum [13] and they represent only 
approximate solutions of the classical equations of motion. This approximation becomes 
better and better as the ratio p/(La) decreases. On the lattice, this problem would lead 
to metastability even in the case of an ideally perfect discretization. However the above- 
defined artificial 1-instanton lattice configurations are enough for our purposes since they 
will be used only as test configurations for cooling methods and not to extract physical 
information. 



In Fig. we show the distribution of the angles 6 z (x) and 6\(x, 1) for an artificial 
1-instanton with size p = 10a located in the middle of a 60 2 lattice. These angles are 
defined as 



|A„(x)-Ar(z)ll=2sin 



0\(x,p) 



\z(x) - z new (x) 



= 2 sin 



x 



(17) 



It is clear that, while the standard cooling acts indistinctly on each lattice site, the 
Pisa cooling and the new cooling deform 1-instanton lattice configurations starting from 
different regions of the lattice. The new cooling will act first on the region around the 
center of the instanton, the Pisa cooling on the border regions. Therefore it is not obvious 
whether under the three coolings the configurations will look the same sweep after sweep. 
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Figure 2: g S L /(2n) (left) and Q L (right) versus p/a during a protracted cooling iteration 
with the three cooling methods. The starting configuration is a 1-instanton with size 
p = 10a placed in the middle of a 60 2 lattice. 



We did therefore the following check: we performed a long (0(1000) steps) iteration of the 
standard and Pisa coolings^ on an artificial 1-instanton configuration with size p = 10a 
and after each sweep we extracted the topological charge Q L , the action S L and the size 
p/a. This size was determined by a very local procedure, namely by looking for the 
maximum of the lattice action density s^ ax and using the relation s^ ax = 2a 2 /p 2 , valid 
for a 1-instanton configuration. Then, starting from the same configuration, we applied 
the new cooling for 8 as small as 5 x 10~ 4 and for each configuration obtained during 
the cooling iteration, we determined again Q L , S L and p/a. Finally, we put on a plot 
Q L and S L as functions of p/a. The result is shown in Fig. |^: the values of Q L and S L 
during the three cooling methods lie on the same curve. The same plot was obtained by 
using a different procedure to determine the size of the cooled configurations, namely a 
global fit of the lattice action density to the continuum expression for the action density 
of a 1-instanton configuration, Eq. QUI). We interpret these results as an indication that 
the three coolings deformed in the same way the starting configuration. If this conjecture 
is correct, it should be possible to put into correspondence the 8 parameter of the new 
cooling with the number of steps of the standard or Pisa coolings. This will be done in 
Section []. 



5 Cooling I— A pairs 



An I-A configuration in the continuum 2d CP^ 1 models can be taken in the form given 

w(x) 



in Refs. 15, tt 



z\x) 



\W(X) 



= iz{x) ■ d^x) 



1 For the Pisa cooling we always use i5pi sa = 0.2, following Ref. p3 
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log E log 



Figure 3: G(d) (left) and G(p) (right) versus the internal energy E during a pro- 
tracted cooling iteration with the three cooling methods in a log-log scale. The start- 
ing configuration is an I-A pair with centers = (15a, 30a), = (45a, 30a) (hence 
d = — c^l = 30a) and p = 5a on a 60 2 lattice. The arrows show the direction of the 
evolution under cooling. 



with 

w 1 (x 1 ,x 2 ) = (xi + ix 2 — ai)(xi — ix 2 — h) , 

w 2 (x 1 ,x 2 ) = (xi + ix 2 - a 2 )(x 1 - ix 2 - b 2 ) , 

Wi = 0, i = 3,...,N. (19) 

In this expression, the complex numbers a\ t2 and b\ t2 are related to the position of the 
center and the size of the instanton and anti-instanton through the following relations: 

a, = 4 7) + ic? - pV , a 2 = cS 7) + ic? + , 

61 = ci A) + ic^ - pW , b 2 = c[ A) + 4 A) + pW , (20) 

where c[ f ' A ^ is the i-th coordinate of the center of the instanton or anti-instanton and 
p(i,A) are their sizes. As "artificial" I-A lattice configuration we took the discretization of 
the above continuum configuration on a 2d finite volume lattice with periodic boundary 
conditions. For simplicity, we used = p^ = p. 

Starting from several I-A pairs with different values of the distance between the centers 
and and of the size p, we performed long cooling sequences with the three methods, 
namely 0(1000) iterations for the standard and Pisa coolings and 5 as small as 5 x 10~ 4 
for the new cooling. In order to have an indication of how the cooling process works, we 
examined by eye the distribution of the topological charge and action densities during the 
different cooling iterations. The evolution of these densities looked the same, no matter 
which cooling procedure was adopted. Specifically, the distribution of the action density 
shows at the beginning two equal instanton bumps, which for d ^> p (d is the distance 
between the two centers and c^) are well-separated; then, as the cooling goes on, 
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these bumps lower in height and merge together, up to complete annihilation. In order to 
make quantitative the statement that the three coolings behave essentially in the same way, 
we studied the shell correlation function of the topological charge density. This function 
is defined as G(r) = l/N r ^2 xy q L (%)q L (y) where the sum is extended over all pairs of 
sites x, y which satisfy r < \x — y\ < r + Ar and N r is the number of those pairs for a 
given value of r (we have chosen Ar = 0.6 a). We determined G(r) on every configuration 
resulting after each cooling iteration^. We compared the results which presented the same 
value for the internal energy E which is defined by 

E = ^2 Yl + • z ^ M*) + z ^ ■ z ( x + M x ) - 2 ) > ( 21 ) 

x n=l,2 

and is proportional to the lattice action. In Fig. ^ we show the behavior of the shell 
correlation versus E during the three (protracted) coolings for the case of r equal to the 
distance d between the centers of the instanton and of the anti-instanton (left) and to the 
instanton (or anti-instanton) size p (right). The three curves fall on top of each other, 
apart from a small deviation between the Pisa cooling and the other two at the very first 
(1-2) cooling steps, due to (unphysical) finite size effects. Since the shell correlation of 
the topological charge density is related to the size of the instantons, the conclusion is 
that the three cooling methods perform equivalently in deforming the size of both the 
instanton and the anti-instanton in the I-A pair and in modifying the distance between 
the two peaks. 



6 Cooling on the equilibrium ensemble 

Using artificial configurations with non-trivial topology provides interesting suggestions 
about a cooling procedure, but in real life the cooling must be applied on thermal equilib- 
rium configurations. Among them there are configurations which contain several instanton 
and anti-instanton bumps, merged in a sea of quantum fluctuations. These configurations 
can represent as well good test-fields for the different cooling procedures. Their evolu- 
tion under cooling can be visualized by the distribution of the action density or of the 
topological charge density. In principle we could expect that the three cooling methods 
described in Section [3] make the starting thermalized configuration evolve along different 
directions. In particular, for a not too protracted cooling, one expects that all the three 
methods have been able to erase the quantum fluctuations, but the number, type and 
location of the surviving topological bumps in the cooled configuration can be rather dif- 
ferent. In order to check this expectation, it is necessary to find a criterion to make a 
correspondence between the number of iterations in the standard or Pisa coolings with a 
5 parameter for the new cooling. We defined an effective temperature for configurations 
obtained during the cooling iteration in such a way that the comparison could be made 

2 From reflection positivity it follows that the two-point correlation function of the topological charge 
density, and consequently also (G(r)), is negative at distances r > 0. However reflection positivity is lost 
during cooling, since cooling affects the quantum structure of the theory leaving only the semiclassical 
background intact. After cooling one in general obtains positive values for (G(r)), apart from those large 
separations which the cooling procedure has not yet affected. 
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Figure 4: Comparison of the determinations of \ L after equivalent amounts of new and 
standard coolings for the CP 3 model at (3 = 1.05 on a 76 2 lattice. The new cooling results 
are reported versus equivalent iterations of standard cooling steps following Table 1. Full 
(open) symbols indicate the new (standard) cooling. 



only between configurations at the same temperature. The most natural "thermometer" 
is the internal energy E, Eq. (^1|), since this is the quantity which is minimized during the 
cooling procedure. 

We put into practice the above considerations in the case of the CP 3 model. We 
generated by Monte Carlo a sample of equilibrium configurations on a 76 2 lattice at (3 = 
1.05. As in Ref. |14j], the simulation algorithm is, for every updating step, a mixture of 
4 microcanonical updates and 1 over heat-bath [17|. We measured the expectation value 



of E on the cooled ensembles obtained by the three cooling methods after several number 
of iterations (standard and Pisa coolings) and for several values of 5 (new cooling). By 
comparing the results and imposing that E be the same, we obtained a correspondence 
table between the number of cooling iterations for the standard and Pisa coolings with 
a value of the 5 parameter in the new cooling. We found for instance that for the CP 3 
model at the above values of (3 and lattice size, 30 iterations of the standard cooling 
correspond (in the sense of the average value of E) to 5 = 0.007 for the new cooling and 
to approximately 33 iterations of the Pisa cooling. 

In general the matching between the number of standard and Pisa cooling iterations 
is not as good as the one obtained between 5 in the new cooling and a fixed number of 
iterations in the standard (or Pisa) cooling. The reason is that in the second case one 
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Figure 5: The same comparison as in Fig. |] for the magnetic susceptibility and the cor- 
relation length as defined in Eq. (f^) . Data points have been horizontally shifted for the 
sake of readability. Full (open) symbols indicate the new (standard) cooling. 



can tune a continuous parameter while in the first only integer jumps of the parameter 
controlling the amount of cooling are allowed. For this reason we have mainly investigated 
the relation existing between the standard and new coolings. In Table 1 we show the 
values of S (new cooling) found in correspondence to several numbers of iterations of the 
standard cooling for the CP 3 model at (3 = 1.05 on a 76 2 lattice. Notice that a different 
correspondence table may be found for other values of iV and f3. 

Using the correspondence table we can directly compare the different coolings. As a 
first step, we have computed the average values of some physical quantities on samples ob- 
tained after applying equivalent amounts of cooling on a set of equilibrium configurations. 
In Fig. £| we show the results for the lattice topological susceptibility with the standard 
and new coolings. By using the relation given in Table 1, we report both determinations 
in terms of the number of standard cooling iterations. It clearly appears that, regard- 
ing to the determination of the topological susceptibility, the two cooling techniques are 
completely equivalent. 

Analogously in Fig. [5] we compare the determinations of the magnetic susceptibility 
Xm and of the correlation length £ G defined by 



Xr. 



a 



-l^(TrP(x)P(y)) con 
G P (p,0) 



x,y 
1 



4sin 2 (7r/L) 



g p (o,i; 



(22) 



where 



G P {k) = ^]T(TrP(x)P(?/)> conn exp 



x.y 



4<*-»)-* 



(23) 



is the lattice Fourier transform of the correlation of two local gauge-invariant composite 
operators Pjj(x) which have been defined in Eq. (I7|) (the subscript "conn" means connected 
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Figure 6: Two-point shell correlation of the lattice topological charge density operator (r 
is the distance between two lattice sites) after cooling for the CP 3 model at f3 = 1.05 on 
a 76 2 lattice. The data have been slightly horizontally shifted. 

Green function). The results for Xm and ^g/ - from Fig. |5] indicate that in this case the new 
cooling and the standard cooling show a good agreement, apart from a small discrepancy 
at large numbers of iterations. 

We have also studied the shell correlation of two lattice topological charge density 
operators, (G(r)), for r/a going from to 20, after equivalent amounts of the three coolings. 
In Fig. |] we compare the standard cooling (30 iterations) to the new cooling (5 = 0.007) 
and the Pisa cooling (33 iterations, 5p isa = 0.2)0. The results are in good agreement, 
although a small deviation can be observed at small distances. The suggestion of this 
outcome is that the three coolings affect the instanton size distribution in the small distance 
region in a slightly different way; the effect is however modest. 

Next we compared the three coolings on a subset of the above-generated thermal 
equilibrium ensemble for the CP 3 model. To this aim, we have chosen from the thermal 
ensemble several configurations showing non-zero topological charge after 30 iterations of 
the standard cooling and compared their action density and topological charge density 
distributions by eye to those of the same configurations cooled by 33 iterations of the Pisa 
cooling and by the new cooling with 5 = 0.007. For all the configurations considered, we 
observed that the bumps corresponding to instantons and anti-instantons had the same 

3 See footnote Bin Section about the positivity of (G(r)) after cooling. 



12 




Figure 7: Action density for a configuration obtained by cooling a thermal equilibrium 
configuration of CP 3 at (3 — 1.05 on a 76 2 lattice by the new cooling with 5 = 0.007 
(left) and by 30 iterations of the standard cooling (right). The two amounts of cooling are 
equivalent in the sense explained in the text. 

shape, number and location (see Fig. |7| for an example). 

These considerations provide evidence that the three cooling methods defined in Sec- 
tion |3| behave essentially in an equivalent way, not only on average, but also configuration 
by configuration, and that they can be related by a simple correspondence between number 
of iterations (standard and Pisa coolings) and 6 (new cooling). 

We close this Section with some considerations about the usual picture of cooling as 
a diffusion process. As we pointed out in Section ^ the parameter S in the new cooling 
behaves, in the continuum limit, like a physical quantity of dimension two. Therefore, 
indicating with r the physical scale up to which cooling affects the quantum fluctuations 
and following Ref. |J, we infer that r$ oc <5 -1 / 2 in the continuum limit. On the other 
hand, the standard cooling is usually believed to act as a diffusion process, so that, if n is 
the number of iterations, it should affect fluctuations up to a scale r oc n 1 / 2 . Using the 
relation found between n and 5, we can check both these predictions: indeed we expect 
that 5 oc 1/n. In Fig. || we have plotted the function 5(n) given in Table 1, together 
with a best fit to a function proportional to 1/n. For the lattice parameters used in the 
simulation we see that 8(n) behaves as expected for small values of n, while the relation 
5 oc 1/n breaks down for n > 30 — 40, indicating that the picture of the standard cooling 
as a diffusion process may work well only for moderately small values of n. 



7 The topological susceptibility 



In this Section we apply the new cooling method to the determination of the physical 
topological susceptibility and compare the results to those from an alternative approach, 
the field theoretical method [Ti| combined with smearing [HI]. As for the comparison 
of the new cooling with the standard and Pisa coolings, we have already shown in the 
previous Section that they give equivalent results for \ L ■ 
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Figure 8: Relation between the new cooling parameter 5 and the number of standard 
cooling iterations n, at /3 — 1.05 on a lattice 76 2 in the CP 3 model, together with a best 
fit to a function proportional to 1/n. 

The lattice discretization of the topological charge density and of the topological sus- 
ceptibility has already been discussed in Section 0. In the relation existing between the 
lattice topological susceptibility \ L an d the continuum one, Eq. ([[(]), M(j3) indicates the 
mixings of \ L to operators with the same quantum numbers. More precisely one can write 

M(P) = A(P) a 2 (T) np + P(P)(I) + 0(a 4 ) . (24) 

The first and second terms in the r.h.s. are the mixings with the trace of the energy- 
momentum tensor and with the unit operator respectively ("np" means the purely non- 
perturbative part of (T)). 

The mixing coefficients A(P) and P(P) as well as the multiplicative renormalization 
Z(P) can be calculated in perturbation theory. The perturbative series for A{(3) starts 
from the order 1//3 3 and arguments can be given [|ll|] which justify that A(j3)a 2 {T) np is 



safely negligible in the scaling window of the simulation. In Ref. |2t| the perturbative 
series for Z(P) has been calculated up to the order 1//3 2 and the perturbative series of 
P(P) at the first two non-zero orders 1//3 4 and 1//5 5 . On the other hand there are no 
available perturbative estimates of the 0(a 4 ) terms in Eq.(24). 



However, a more powerful, purely numerical technique can be used to get a non- 
perturbative determination of Z(P) and of the whole additive renormalization M((3). This 
technique is the so-called "heating method" [ZT|. 



The idea of the heating method is to determine the average values of topological quan- 
tities on samples of configurations obtained by thermalizing the short-range fluctuations, 
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Figure 9: x^ooi as a function of the 5 parameter of the new cooling for CP 9 at (3 = 0.70 on 
a 24 2 lattice. 

which are responsible for the renormalizations, on configurations of well defined topolog- 
ical background. For instance if one applies several thermalization steps at a given value 
of (3 on a configuration containing one discretized instanton of charge Q = 1 and measures 
the average value of Q L , then the renormalization constant Z(f3) can be determined as 
Z(j3) = (Q L )/Q (here (•) indicates an average in the background of a fixed charge Q). 
Analogously, by thermalizing a trivial configuration (for example a configuration where 
for all lattice sites x, the respective fields are z(x) = (1,0,0,...) and A M (x) = 1) at a 
given (3 and measuring the average value of \ L i one gets a determination of M(f3). In 
the first case one obtains a non-perturbative estimate of Z(/3) as long as the MS value 
of the background topological charge Q is chosen (see Section ^|). In the second case the 
method amounts to impose the requirement that the physical susceptibility x vanishes in 
the absence of instantons. 

The heating method has been extensively applied in the following. Once all the renor- 



malization constants are known, the extraction of \ from Eq. (10) is possible and the 
result is called the field theoretical determination. 

Since the continuum \ is extracted from the lattice \ L by subtracting the renormaliza- 
tion effects, it happens that, if M(/3) is a large part of the whole lattice signal and Z(/3) is 
small, the physical signal for \ is extracted with very large errors. This is the case when 
one uses the lattice discretization of Eq. (H). However, since both Z(j3) and M(/3) depend 
on the discretization q L (x), one can exploit the arbitrariness in the lattice definition and 
use an improved operator. To this aim, following the idea of Ref. [pT9|| , already used for the 
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determination of x i n Yang-Mills theories [22], we have used a smeared topological charge 
density operator, which is built from the original operator q L (x), defined in Eq. @, by 
replacing the fields z(x) and A M (x) with 



z~(x) 



[1 - c)z(x) + ~ ^2 {[ z ( x - V) K( x - V) + ^ + 

M=l,2 



A— (x) = A/a 



;i- C )A^) +C || gg;;g+g M ], (25) 

where J\f z , A/a are normalization constants which allow that ||z smear (x)|| = ||A^ mear (a;)|| = 1 
and c has been chosen to be equal to 0.65 The smearing can be iterated at will by 
defining the n-th level of smeared fields from the (n — l)-th level in an analogous fashion 
as shown in Eq. fl25|). For each level of smearing a relation like Eq. ( |TUD holds although 
Z(/3) and M{j3) change. Using smeared operators these renormalization constants get 
closer to 1 and 0, respectively, as the smearing level increases, thus allowing a much more 
precise determination of x- 

We performed a numerical simulation for the CP 9 model at several values of (3 and 
L. We used the same updating algorithm of the previous Section. For most simulations 
we collected 100K equilibrium configurations after discarding 10K configurations to al- 
low thermalization. The successive equilibrium configurations were decorrelated by 10 
updating steps. By averaging on the thermal equilibrium ensemble we determined the 
unsmeared lattice topological susceptibility x L an d the smeared ones, Xsmean f° r 1 to 10 
smearing levels. These results, together with Z(j3) and M((3) coming from the heating 
method, allow to extract a?x from Eq. ([TC|) . 

Every 50 updating steps we have also applied the new cooling method at several values 
of the 5 parameter ranging from 0.3 to 5 x 10~ 4 . Since cooling eliminates the short- 
distance fluctuations, the value of x L after cooling, hereafter called Xcooi> can be written 
as in Eq. (|1(]) with Z m 1 and M ~ if the cooling has been protracted enough. Hence it 
should provide directly a?x- Consistency requires that the cooling and the field theoretical 
methods provide the same result. 

The summary of the performed simulations is presented in Table 2. In Fig. |9| we show 
the results for Xmoi fr° m the new cooling at several values for the S parameter at (3 = 0.70 
on a 24 2 lattice. We interpret the behavior under the new cooling in the following way: 
as 5 decreases, the cooling erases an even larger amount of quantum fluctuations thus 
bringing xfooi to a maximum. Afterwards, for too small values of 5, the cooling begins to 
affect also the topological fluctuations and part of the topological signal gets lost. These 
considerations suggest that a convenient choice of the 5 parameter may be around the 
value for which Xcooi reaches the maximum. At this value of 5 we then assume that Z«l 
and M w 0. The choice of a smaller value for S would not be dramatically dangerous, but 
would push the scaling region towards larger values of f3. 

4 This choice has been based on the fact that i) it is convenient to take c as large as possible and ii) 
for c > 2/3 the smearing procedure behaves in a radically different way, leading, when iterated, to a 
completely disordered field configuration instead of a smoother one. A similar behavior has been observed 
in Yang-Mills theories p3| . 
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Figure 10: Scaling of the topological susceptibility of CP 9 from the field theoretical with 
smearing (left) and from the new cooling (right). Data points have been slightly shifted 
along the (3 axis for the sake of readability. 



In Table 3 we summarize the results for \ L an d for £q \ obtained for CP 9 by the field 
theoretical method combined with smearing (0, 5 and 10 levels). In Table 4 the results for 
X^ 00 i and x obtained by the new cooling method for 5 around the peak value (see Fig. |9| 
for p = 0.70) are shown. Fig. displays the scaling of £qX for the cases of 0, 5 and 10 
smearing levels and the field theoretical method (left) and with the new cooling method 
for 5=0.0100, 0.0050, 0.0025 (right). From this figure we see that there is practically no 
dependence on the smearing level, although error bars strongly decrease with the number 
of smearing levels. As for the cooling method, we see that there is a satisfying consistence 
among the results for the different values of S and also between these results and those 
from the field theoretical method. 



8 Summary and conclusions 

The study of the topological properties of the vacuum of a field theory simulated on a lat- 
tice is made difficult by the presence of quantum fluctuations which hinder the extraction 
of the relevant physical information. Among other methods, cooling has been proposed as 
a technique to wash out such fluctuations and reveal the topological background of any 
single field configuration. In the present paper we have made a comparison among three 
different cooling methods using as test-field the 2d CP^ -1 model defined in Section 0. 

The three cooling methods under study have been defined in Section |3|. They are: the 
standard cooling (a local minimization of the Euclidean lattice action), the Pisa cooling 
(the same but with local minimizations constrained to be smaller than a given bound) and 
a new cooling (recently introduced in Ref. ||) where the local minimizations are accepted 
only if they are larger than a given bound. 

Firstly, we have studied the performance of the three coolings on classical non-equi- 
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librium configurations representing 1-instanton solutions and instanton-anti-instanton 
pairs. We have placed one such object on the lattice and have performed a series of 
cooling iterations in order to measure several physical (topological and non-topological) 
observables. In all situations the three coolings have yielded the same result — see Figs. |2] 
and || 

After these results, we have assumed that the three cooling methods are equivalent and 
that a correspondence can be established between them in the sense that a given number 
of iterations of the Pisa or the standard coolings corresponds to a precise value for the 
bound in the new cooling if the value of the energy E, Eq. fl2lD , is the same after applying 
the coolings. This correspondence has been used in the subsequent investigation and, for 
the CP 3 model at /3 — 1.05 on a 76 2 lattice, it is shown in Table 1. 

We have studied the performance of the cooling methods on a set of equilibrium config- 
urations obtained after a Monte Carlo simulation. On these thermalized configurations we 
have first extracted several physical quantities after applying cooling: the lattice topolog- 
ical susceptibility is the same when equivalent amounts of cooling, in the sense of Table 1, 
are applied — see Fig. [|. The same happens for the magnetic susceptibility and for the 
correlation length as far as the cooling iteration is not protracted too much — see Fig. [5]. 
The same agreement is seen for the shell correlation function (G(r)) if r/a is large enough 
— see Fig. |[ Then we have peered the action density and the topological charge density 
distributions obtained with the application of equivalent amounts of cooling, in the sense 
of Table 1. In all cases we have again obtained completely analogous distributions — see 
Fig. [7] for an example. These results altogether strongly suggest that the three cooling 
methods, although different in the procedure, behave equivalently. 

By using the correspondence given in Table 1 we have tested the picture of cooling as 
a diffusion process. For the lattice parameters used in this analysis, this picture works 
well only for a moderately small number of iterations (n < 40). 

Finally we have compared the results obtained for the topological susceptibility in the 
CP 9 model by the new cooling method with those extracted from a well-tested method: 
the so-called field theoretical method improved with smearing. In Tables 3 and 4 we give 
the results for the quantity £q X- They agree fairly well among them, see Fig. [It], and also 
with the large N estimate |24[] which provides £qX — 153 up to order 0(1/N 2 ). 

We used the standard action both to generate thermal equilibrium configurations in 
Monte Carlo simulations and during the cooling procedure. We expect that the above 
results concerning the equivalence between the different cooling techniques should not 
depend strongly on the action used during cooling. However the comparison of the three 
cooling techniques by using different lattice actions is worth to be pursued in a future 
work. 
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TABLES 



Table 1: Relation between the new cooling parameter 8 and the number of standard 
cooling iterations n, at f3 — 1.05 on a 76 2 lattice in the CP 3 model. The errors given in 
the second column derive from the uncertainty in the matching procedure. 

n (standard cooling) 5 (new cooling) 

3 0.068240(50) 

4 0.048500(50) 

5 0.038000(40) 
7 0.026800(30) 
10 0.018840(40) 
15 0.012850(40) 
20 0.009925(40) 
30 0.007000(30) 
40 0.005492(20) 
50 0.004602(20) 
60 0.003990(30) 
70 0.003530(30) 
100 0.002720(20) 
150 0.002038(10) 

200 0.001662(10) 



Table 2: Summary of the simulations for the CP 3 and CP 9 models. 



N 


p 


L 


stat 


E 


id a 


4 


1.05 


76 


100k 


0.530638(12) 


7.688(29) 


10 


0.70 


24 


100k 


0.78429(3) 


2.316(3) 


10 


0.75 


30 


10k 


0.72013(8) 


3.297(11) 


10 


0.75 


32 


100k 


0.720146(25) 


3.2848(24) 


10 


0.75 


60 


10k 


0.72024(5) 


3.248(25) 


10 


0.80 


48 


100k 


0.667000(15) 


4.608(4) 


10 


0.85 


64 


100k 


0.622271(10) 


6.419(7) 


10 


0.85 


80 


10k 


0.622357(29) 


6.371(29) 


10 


0.90 


90 


100k 


0.583830(7) 


8.836(10) 
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Table 3: Summary of the results for the topological susceptibility of the CP 9 model by 
the field theoretical method at various smearing levels. 





L 


level 


X L x 10 4 


M{(5) x W 


Z{0) 


a 2 xx 10 4 


CgX* 10 4 









0.992(5) 


33(3) 


0.160(10) 


26(5) 


139(25) 


0.70 


24 


5 


16.72(8) 


78(15) 


0.810(20) 


24.3(1.6) 


130(9) 






10 


21.40(11) 


45(15) 


0.907(15) 


25.5(1.2) 


137(7) 









0.813(12) 


23.0(1.0) 


0.197(10) 


15.0(2.1) 


163(24) 


0.75 


30 


5 


10.41(16) 


16.0(2.5) 


0.837(15) 


14.6(8) 


159(10) 






10 


12.60(19) 


5.5(1.0) 


0.915(15) 


15.0(7) 


163(9) 









0.809(4) 


23.0(1.0) 


0.197(10) 


15.0(2.1) 


161(20) 


0.75 


32 


5 


10.43(5) 


16.0(2.5) 


0.837(15) 


14.6(8) 


158(7) 






10 


12.67(6) 


5.5(1.0) 


0.915(15) 


15.0(7) 


163(6) 









0.810(12) 


23.0(1.0) 


0.197(10) 


15.0(2.1) 


158(24) 


0.75 


60 


5 


10.51(15) 


16.0(2.5) 


0.837(15) 


14.8(8) 


156(11) 






10 


12.82(19) 


5.5(1.0) 


0.915(15) 


15.2(7) 


161(10) 









0.5981(20) 


16.0(1.5) 


0.240(10) 


7.6(9) 


161(20) 


0.80 


48 


5 


5.87(3) 


2.0(6) 


0.876(8) 


7.63(19) 


162(4) 






10 


6.87(4) 


0.40(15) 


0.933(7) 


7.88(17) 


167(4) 









0.4208(28) 


13.0(7) 


0.270(6) 


4.0(3) 


164(13) 


0.85 


64 


5 


3.180(27) 


0.40(20) 


0.892(7) 


3.99(10) 


164(4) 






10 


3.61(3) 


0.05(3) 


0.945(7) 


4.05(9) 


167(4) 









0.423(8) 


13.0(7) 


0.270(6) 


4.0(4) 


163(17) 


0.85 


80 


5 


3.19(8) 


0.40(20) 


0.892(7) 


4.00(17) 


162(8) 






10 


3.63(9) 


0.05(3) 


0.945(7) 


4.06(17) 


165(8) 









0.2975(15) 


9.5(5) 


0.284(7) 


2.51(20) 


196(16) 


0.90 


90 


5 


1.74(3) 


0.10(3) 


0.902(8) 


2.14(7) 


167(6) 






10 


1.94(3) 


0.010(3) 


0.945(9) 


2.18(8)) 


170(6) 
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Table 4: Summary of the results for the topological susceptibility of the CP 9 model from 
the new cooling method at various 5 values. 





L 


5 


Xcooi x 10 4 


^XXlO 4 






0.0500 


24.1(3) 


129.0(2.1) 






0.0250 


26.7(4) 


143.3(2.4) 


0.70 


24 


0.0100 


28.6(4) 


153.3(2.6) 






0.0050 


28.4(4) 


152.2(2.7) 






0.0025 


26.5(4) 


142.1(2.6) 






0.0500 


13.5(5) 


146(6) 






0.0250 


14.8(5) 


161(6) 


0.75 


30 


0.0100 


15.9(5) 


173(7) 






0.0050 


16.2(5) 


177(7) 






0.0025 


16.1(5) 


175(7) 






0.0500 
0.0250 


13.70(20) 
15.08(22) 


147.8(2.4) 
162.7(2.6) 


0.75 


32 


0.0100 


16.19(24) 


174.7(2.8) 






0.0050 


16.58(24) 


178.9(2.9) 






0.0025 


16.51(25) 


178(3) 






0.0500 


13.5(4) 


143(7) 






0.0250 


14.8(5) 


156(7) 


0.75 


60 


0.0100 


15.9(5) 


167(8) 






0.0050 


16.4(5) 


173(8) 






0.0025 


16.6(5) 


175(8) 






0.0500 


7.12(9) 


151.1(2.1) 






0.0250 


7.73(9) 


164.1(2.3) 


0.80 


48 


0.0100 


8.22(10) 


174.6(2.5) 






0.0050 


8.44(10) 


179.2(2.5) 






0.0025 


8.57(10) 


181.9(2.5) 






0.0500 


3.68(4) 


151.8(1.9) 






0.0250 


3.94(4) 


162.5(2.0) 


0.85 


64 


0.0100 


4.14(4) 


170.7(2.0) 






0.0050 


4.23(4) 


174.4(2.1) 






0.0025 


4.29(4) 


176.8(2.1) 






0.0500 


3.63(12) 


147(6) 






0.0250 


3.89(13) 


158(7) 


0.85 


80 


0.0100 


4.09(14) 


166(7) 






0.0050 


4.17(14) 


169(7) 






0.0025 


4.23(14) 


172(7) 






0.0500 


1.97(3) 


153(3) 






0.0250 


2.08(4) 


162(3) 


0.90 


90 


0.0100 


2.17(4) 


169(3) 






0.0050 


2.21(4) 


172(3) 






0.0025 


2.23(4) 


174(3) 



24 



